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Abstract. We report on four observations of the black hole candidate GX 339-4 collected in the spectral range 0.4-140 keV 
with the BeppoSAX satellite during 1997. All our spectra are typical of a black hole candidate in its low-hard state, i.e. are 
dominated by the hard tail of the emission from the corona. Although the direct emission from the accretion disk in the soft 
X-ray range is weak, spectra are well described only with models that take into consideration the Compton reflection of hard 
X-rays from the corona onto the surface of an accretion disk of non negligible extent. 
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1. Introduction 

GX 339-4, discovered in 197 1 by Markert et al. <1973alfT973bl . 
was proposed to be a black hole candidate by Samimi et al. 
( 1979 1, because of its spectral features and its short term vari- 
ability similar to sources like Cyg X-l. 

GX 339-4 was observed in five different spectral states, 
each one with well defined features (Tanaka & Lewin 1995 1. 
Most observations detected the source in the hard-low state, 
and often with such a low count rate to properly call it an off 
state; less often the source was detected in the high-soft state, 
only twice in the very high state (before 2001) and twice in 
the intermediate state (for a recent re-consideration of the na- 
ture of the intermediate state see Homan et al. 12001b . In the 
first half of 1999, the source entered an extended off-state, 
and it remained undetectable by the All-Sky Monitor on board 
RossiXTE (see Kong et al. 2000 Corbel et al. 120031 . During 
2002 GX 339-4 started a new outburst (Smith et al. l2"0"02T > with 
more state transitions (Belloni et al. 2003 1 The light curve from 
the RXTE/ASM is shown in Fig.Q] 

In its low state GX 339-4 shows power law spectra, with 
photon index r=1.62-2.30 (Dolan et al. 119871 Bouchet et al. 
\TW5\ Ueda et al. [T9^4l Mendez & van der Klis flW7l Zdziarski 
et al. fT998l Belloni et al. fT999l Wilms et al. [1999). Higher 
values for F were detected in spectra observed up to 200 keV 
(F = 2.0 + 0.3; Dolan et al. [19871 and to 1.3 MeV (Bouchet et 
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al.[T993 1. Most observations in the low state indicate the pres- 
ence of the fluorescence Fe line (Ueda et al. 1 19941 Zdziarski et 
al l 19981 Dolan et al. U987> . with an equivalent width of about 
50-100 eV. 

Several authors (Zdziarski et al. IT9981 Wilms et al. IT9991 
Revnivtsev et al. 120011 Nowak et al. 2002 Wardziriski et al. 
2002 1 give indication of a reflection component in the X spec- 
tra, described according to the model of Magdziarz & Zdziarski 
( 11995b and provide evidence for the Fe abundance being greater 
by a factor 2-3 relative to solar abundance. 

GX 339-4 high states have also well defined spectral fea- 
tures. In this state the spectrum shows a soft component, due 
to the thermal emission by an accretion disk (Makishima et al. 
1986 Belloni et al. 1999 1, described by a multicolour black- 
body model (Makishima et al. 1986 1. Inner temperatures of the 
disk were estimated to be 0.77 keV (Makishima et al. 1 1986 1 
and 0.63-0.72 keV (Belloni et al. [T999l . The high energy part 
of the observed spectra has a powerlaw behaviour, with a pho- 
ton index T of 2.0 (Makishima et al. [T9861 and of 2.12-2.57 
(Belloni et al. fT999l 

In the high state the fluorescence Fe line at 6.4 keV is 
present, with an equivalent width of 319-350 eV. There is also 
evidence for an absorption edge, although its central energy is 
uncertain: Makishima et al. (1986 1 report a value of 8.8 keV, 
while Belloni et al. fl999l one of 6.69 keV. 

The source is very active also in the radio band, with radio 
flux strictly correlated to the X-ray intensity (Corbel et al. 2003 
Markoffet al. l2003l . 
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Fig. 1. Historical GX 339-4 light curve as observed by the All Sky Monitor onboard Rossi XTE. The dashed vertical lines in the 
main panel and in the inset indicate the epoch of our BeppoSAX observations. 



In this paper we report results from four observations of 
GX 339-4, taken during 1997 with the BeppoSAX satellite, 
when the source was in a low-hard state. Using all four narrow 
field instruments available in the satellite scientific payload, we 
obtained and analysed data in the 0.4-140.0 keV band. In Sect. 
2 we present the four observations with a summary of the re- 
duction procedure; in Sect. 3 we describe the data analysis and 
in Sect. 4 we discuss our results. 

2. Observations and data reduction 



We observed GX 339-4 with BeppoSAX (Boella et al. IT9971 
and references therein) at four epochs during 1997, with the 
primary aim to study the spectral variability on the 1 -month 
and 6-month time scale. A journal with basic data is reported 
in Table [0 while the state of the source at the epoch of our 
observations with respect to the historical light curve is visible 
in Fig.[0 Preliminary results on the first two observations have 
been previously reported by Chiappetti et al. ( 119981 . 

From the figure it is apparent that we observed the source 
at only slightly different intensity levels during a low state. 

We detected GX 339-4 with all the four Narrow Field 
Instruments on board BeppoSAX, obtaining a simultaneous 
coverage with an excellent signal to noise from 0.4 to 140 keV. 
MECS (Medium Energy Concentrator Spectrometer) observa- 
tions used all three MECS units for the first two observations, 
and only two units for the last two, due to a failure of unit 
MECS1 in late spring 1997. 



Data reduction for MECS, HPGSPC (High Pressure 
Gas Scintillation Proportional Counter) and PDS (Phoswich 
Detector System) was performed using the XAS software 
(Chiappetti & Dal Fiume [T997> . while for LECS (Low Energy 
Concentrator Spectrometer) data we used the cleaned event 
files supplied by the ASI Science Data Centre (ASDC). In gen- 
eral the reduction followed the standard prescriptions described 
e.g. in Chiappetti et al. ( 1999 1, with only slight differences and 
peculiarities due to the higher intensity of GX 339-4. A sum- 
mary is given below. The reader not interested in reduction de- 
tails may wish to proceed directly to Sect. 3. 

We examined the light curves of GX 339-4 during our (rel- 
atively short) observations, but, since we found no indication of 
large variability, we proceeded accumulating a single spectrum 
per instrument for each observation. 

LECS gross spectra were extracted from event files in a cir- 
cular area of radius 8 arcmin. Background subtraction was per- 
formed linking each data file to a standard background file (the 
latest available at ASDC 1 in November 1999). ASDC standard 
response matrices were also used. In our analysis we used data 
in the 0.4-4.0 keV energy band. 

A larger extraction radius of 1 3 arcmin was used for MECS . 
Background subtraction was performed using a 1 120 ks back- 
ground file derived from blank field pointings taken during 
the satellite Performance Verification phase. Response matri- 
ces were generated with XAS according to the source location 



calibration files repository: ftp . asdc . asi . it/pub/ sax/ cal/ 
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Table 1. Journal of BeppoSAX observations 



Number 


Start epoch 


LECS 


LECS 


MECS3 


MECS3 


HPGSPC 


HPGSPC 


PDS 


PDS 




(UT) 


exposure 


count rate 


exposure 


count rate 


exposure 


count rate 


exposure 


count rate 


i 


17 Mar 07:08 


4888 s 


6.22±0.04 


10108 s 


6.92+0.03 


4262 s 


21.5±0.3 


4407+4273 s 


18.6+0.1 


2 


06 Apr 10:13 


7302 s 


5.83±0.03 


12286 s 


6.10+0.02 


5628 s 


18.8±0.2 


5631+5626 s 


16.8+0.1 


3 


06 Sep 05:51 


5947 s 


3.96±0.03 


16858 s 


4.72±0.02 


7432 s 


15.7+0.2 


7976+7906 s 


13.9+0.1 


4 


01 Oct 17:36 


5079 s 


9.01±0.04 


25046 s 


9.81±0.02 


11006 s 


27.6+0.1 


10897+11018 s 


22.0±0.1 



a Count rates measured in the energy bands used for the fits, namely: 0.4-4.0 keV (LECS); 4.0-10.0 keV (MECS1); 1.8-10.0 keV (MECS2 

and MECS3); 7.0-70.0 keV (HPGSPC); 20.0-140.0 keV (PDS). 
b For MECS we indicate in the table the count rate for unit MECS3. Unit MECS2 has a similar count rate, while unit MECS 1 (where present, 

i.e. in obs 1 & 2) has a slightly lower count rate due to thicker ion shield, 
c For the PDS instrument we quote the exposure time as the sum of the two instrument halves while the count rate is the net count rate in the 

4 units of which 2 look at the source at any one time for the quoted exposure times 



in the field of view. In our analysis we used data in the 4.0-10.0 
keV energy band for MECS1, in the 1.8-10.0keV energy band 
forMECS2 andMECS3. 

HPGSPC background subtraction was performed using 
data taken during our observations (but not simultaneously) 
with the rocking collimator offset in the negative direction. 
ASDC standard response matrices were used. In our analysis 
we used data in the 7.0-70.0 keV energy band. 

PDS background subtraction used simultaneous data taken 
with the rocking collimators in offset position (Frontera et al. 
119971 . In consideration of the source strength the threshold 
for spurious spike filtering has been 32 cts/s (instead of the 
customary 25 cts/s), and the PSA correction was not applied. 
Response matrices were computed with XAS from the config- 
uration recorded in the spectra. In our analysis we used data in 
the 20.0-140.0 keV energy band. 

Given the brightness of the source, the systematics in the 
calibrations of the instrument play an important role when 
matching spectra taken with different instruments. This is cus- 
tomarily handled using the standard response matrices and 
introducing cross-normalisation constants in the spectral fit. 
Unit MECS3 was used as reference (constant fixed to unity). 
Obtained cross-normalisations for the other instruments are in 
good agreement with values widely reported in the literature 
and recommended by the BeppoSAX team. 

3. Data analysis 

Spectral fits were performed with the software package XSPEC 
v. 10.0. In the remainder of the paper we refer to all additive or 
multiplicative spectral components, as well as their parameters, 
with the mnemonic used in the XSPEC handbook. 

Preliminary fits were made with the model WABS (DISKBB 
+ POWERLAW) , to account for the direct emission from the two 
main structures around a black hole accreting matter from a 
companion in a binary system: the optically thick, geomet- 
rically thin, accretion disk (DISKBB, Makishima et al. [1986 > 
and an optically thin Comptonising corona (POWERLAW). The 
component WABS describes photoelectric absorption due to the 
interstellar matter (Morrison & McCammon 1983 Anders & 
EbiharalHHD. 



GX339-4 17 March 1997 - Model DBB+PO 




1 10 100 

channel energy (keV) 



Fig. 2. Result of the WABS (DISKBB + POWERLAW) fit to the 

spectrum of Mar 17 : unfolded spectrum in top panel, contribu- 
tion to x 1 in bottom panel. 



Values of the x 2 ( see Table |2j indicate that such a sim- 
ple model does not describe satisfactorily the observed spectra, 
and Fig.|2](bottom panel) shows that the main deviations from 
the model are in the high energy tail of the spectra: the model 
underestimates the observed spectra at energies around 30-60 
keV, and overestimates the spectra at higher energies. 

A first possible interpretation of such deviations would be 
that the spectra have a high energy cutoff, which has not been 
taken into account in such a simple model : the underestimate 
at medium energies could be interpreted as the attempt to com- 
pensate the lack of a cutoff in the model. A less important de- 
viation is seen at energies around the fluorescence Fe line. 

Therefore, a new series of fits were performed with the 
model WABS (DISKBB + CUTOFFPL + GAUSSIAN): the com- 
ponent CUTOFFPL was introduced to account for the high en- 
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Table 2. Spectral fits with the WABS (DISKBB + POWERLAW) model. 



Component 


Parameter 


Unit 


17 March 


6 April 


6 September 


1 October 


WABS 

DISKBB 

DISKBB 

POWERLAW 

POWERLAW 


N H 

T, n 

Kdbb 

r 

Kpo 


10 21 cm" 2 
keV 

(km/10 pc) 2 

keV/cm 2 s keV (at 1 keV) 


5 00 

3or!!!! 2 5 

-3979 +2466 

1 69 r - 013 

U - /8 °-0.007 


5.08^ 

35l + !M! 30 

-0.026 

1147 +768 

1 695+ 0013 

n 9C7+O.OO6 
"■^ -0.006 


4 63 +0 - 35 

u -0.30 

360 +0 n ° n 5 ' 

JW -0.040 

501+ 613 

1 666+ 0013 
n r s s +0.004 

U - i88 -0.004 


u -" -0.24 

363+ , , 10 . 
1786.!™; 

1 77^+0.007 
1 - ' '-'-0.007 
n 4^4+0.006 

-0.007 




x 2 




706.78 


791.97 


685.33 


2227.68 




^/D.o.F. 




1.78 


1.97 


1.94 


5.64 



a Errors quoted are at the 90% confidence level. 



GX339-4 17 March 1997 - Model DBB + CUTPL+GAUSSIAN 
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Fig. 3. Result of the WABS (DISKBB + CUTOFFPL + Fig. 4. Result of the WABS (DISKBB + PEXRIV + LA0R) fit 



GAUSSIAN) fit to the spectrum of Mar 17 : unfolded spectrum 
in top panel, contribution to \ 2 m bottom panel. 



to the spectrum of Mar 17 : unfolded spectrum in top panel, 
contribution to^ 2 in bottom panel. 



ergy cutoff in the hard tail of the spectra and the component 
GAUSSIAN for the profile of the Fe fluorescence line. The Fe 
line energy was kept fixed to the value of 6.4 keV (see Table[3j. 
X 1 values improved for all observations, but they were still too 
high to indicate good fits. Fig. [3](bottom panel) shows the pres- 
ence of non negligible modulations, indicating that this im- 
proved model cannot describe satisfactorily the shape of the 
observed spectrum. Fig. [3] (top panel) shows clearly that data 
remain above the model curve in the range about 20-50 keV 
and are definitely below the model in the hard tail. 

The behaviour at energies above 30 keV is suggestive of 
the presence of a reflected component in the spectra. To test this 
scenario, we introduced the additive model PEXRIV (Magdziarz 
& Zdziarski 1995 1, which describes the Compton effect on a 
possibly ionized disk, when the incident spectrum is well de- 
scribed by a power law with a high energy cutoff. In conjunc- 
tion with the PEXRIV component, we performed fits represent- 



ing the Fe fluorescence line with the component LAOR (Laor 
1991 1, which describes the line profile from an accretion disk. 

Fitting results (see Table 0} indicate clearly that the 
Compton reflection scenario gives a good description of our 
data, and Fig. |4] shows how good is the agreement of the as- 
sumed model with the data. 

4. Discussion 

Our results clearly show that the source was in the hard-low 
state during all four observations (see Table |4] for numerical 
results and Fig.|4]for the unfolded spectrum). Count rates and 
luminosities for each observation are all within a factor of two 
from each other (see Tables^andEJ and are typical of GX 339- 
4 in its hard-low state (Zdziarski et al. 119981 Ilovaisky et al. 
IT986I Mendez & van der Klis ll997l Belloni et al. lT999l Wilms 
et al. I1999i . The spectra are also resemblant of black hole can- 
didates in such spectral state, being mainly dominated by a 
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Table 3. Spectral fits with the WABS (DISKBB + CUTOFFPL + GAUSSIAN) model 



5 



Component 


Parameter 


Unit 


17 March 


6 April 


6 September 


1 October 


WABS 

DISKBB 

DISKBB 


N H 

T,„ 

Kdbb 


10 21 cm" 2 
keV 

(km/10 pc) 2 


-X 1 c+0.70 
J -0.29 
f) 7*^7+0.096 
—0.240 

50.24+; 20 , 44 


ri ot+0.26 

3.83_ 025 

664 +0 075 
u.uut „ 064 

81.2+ 4 4 °, 


7 OO+0.18 

f) QQ7+0.043 
" ,7J/ -0.047 

19 4 +33 

7 -2.6 


T 7Q+0.13 

—0. 12 

0.877+°,° 2 2 
72 +6 6 


CUTOFFPL 
CUTOFFPL 
CUTOFFPL 


r 

Ecut 
K CL ]tpi 


keV 

keV/cm 2 s keV (at 1 keV) 


l S() 7+ yl z 

1 77+0.56 

1 ■ ' ' -0.30 


1 47 9 

l 65 +(U6 

" -0.19 


. — — i n a^i 
-0.053 

105.5+; 4 | 

n ss+o-oi 
O.85_ 001 


_ — ^ — l A Al A 

1 262 
74 5 +4 2 

'^• J -3.8 

1 59 +0 12 


GAUSSIAN 
GAUSSIAN 


<T 

Kline 


keV 

phot cnT 2 s _1 (xlO~ 3 ) 


0.84!^ 


65 +0 - 24 

"• UJ -0.19 

1 -^-0.54 


95 +017 

U -^ J -0.17 

2 56 +0 - 76 

^• JU -0.65 


1 07 +00<) 

1 - u/ -0.09 

7 61 +101 




2 

X 




528.90 


546.66 


429.43 


698.87 




^/D.o.F. 




1.37 


1.37 


1.22 


1.79 



a The Fe-line energy En nc in the component GAUSSIAN has been fixed to 6.4 keV 
b Errors quoted are at the 90% confidence level. 



Table 4. Spectral fits with the WABS (DISKBB + PEXRIV + LAOR) model 



Component 


Parameter 


Unit 


17 March 


6 April 


6 September 


1 October 


WABS 

DISKBB 

DISKBB 

PEXRIV 

PEXRIV 

PEXRIV 

PEXRIV 

PEXRIV 


N H 
T,„ 
Kdbb 

r 
p 

Rel.Refl 

F^abund 
Xi 


10 21 cm" 2 
keV 

(km/10 pc) 2 
keV 

Relative to solar abundance 


4 76 +031 

^■'"-0.35 

0.32** 
1847.: 2 «5 
i 70+0.04 

-0.03 
q+461.8 

0.80 + »-^ 
0.37 

(<16.29) 


4 66 +014 
H-u -0.22 

43+ 004 
"^"-0.04 

378.6 +274 » 

1 67 +a04 ' 

1 - u '-0.03 

285.5^°' 6 9 

70 +012 
"■'"-o.io 

4.85!}-| 
4.73xl0~ 2 

(<15.31) 


4 46 +0 06 

^•^"-0.16 
041 +0.04 
"•^ i -0.03 
17Q 7+104.7 
1/5 "-'-40.9 
1 79+O.OI 
L - /z -0.04 

802.7+^-5 
0.88+°;??' 

9 os+0.93 

5.55 x 10~ 12 

(<0.07) 


C QQ+0.14 

J-J -0.13 

49+0.02 

"•^^-0.04 

397.4^ 77 

1 7C+0.01 

L - 1 -0.01 

264.8+21 6 
90+° 04 

"■^"-0.04 

3 50+ 047 

-'• J "-0.53 

9.12xl0~ 6 

(<36.12) 


PEXRIV 
LAOR 


V 

^Pexriv 
KLaor 


keV/cm 2 s keV (at 1 keV) 
photons cnT 2 s -1 (x10~ 4 ) 


O 7Q , I+0.014 
"■" J -fl.OI3 
tt+3.55 
-6.20 


O 7^7+0.013 
-0.006 

5 13 +2 - 84 

-2.03 


195+0-003 

S "31+2.64 
j.oj_ 223 


0416 +0.009 
o.tiu 0017 

1 6 90+ 324 
iu -^"-12.22 




Flux 


erg cnT 2 s~' (xl0~ 9 ) 


5.98 


5.21 


4.34 


7.46 




Luminosity 


ergs" 1 (xlO 37 ) 


2.24 


1.95 


1.63 


2.79 




2 

X' 




391.86 


409.86 


340.87 


446.73 




X 2 /D.o.F. 




1.02 


1.04 


0.99 


1.17 


Cross-normalisation constants with reference to MECS3 = 


unity 










LECS 




0.8145 


0.8065 


0.7921 


0.8139 




MECS1 




0.9768 


0.9758 








MECS2 




0.9723 


0.9731 


0.9806 


0.9743 




HPGSPC 




1.056 


1.044 


1.048 


1.019 




PDS 




0.8116 


0.8492 


0.8096 


0.7834 



a The following parameters have been fixed: Fe-line energy EL aor in the component LAOR to 6.4 keV; the disk inclination angle, present in 

both components PEXRIV and LAOR, to the value cos i =0.5. 
b Errors quoted are at the 90% confidence level. 

c Reported fluxes and luminosities are relative to the 0.4-140 kev energy band. Luminosities assume a distance of 5.6 kpc. 



power law, with a weak soft excess at lower energies (around 
1 keV). The spectral indices of such power laws are also typ- 
ical of hard-low states, and confirm that the spectra are gen- 
erated in an optically thin region, assuming they are due to 
Comptonisation. The large solid angle required by the reflec- 
tion component suggests the presence of an extended accretion 
disk, despite the fact that the direct emission from it accounts 



only for a small fraction of the soft flux (about 3-6% in the 
0.4-10 keV band). 

Very little we can argue from the ionisation parameter: all 
confidence intervals are consistent with zero, but at the same 
time have very high upper limits. Equivalent widths of the Fe 
line (from Tables [5] and 0} are also consistent with previous 
value in litterature (Ueda et al. 1994 Zdziarski et al. 119981 
Dolan et al. U987> . lying in the range 50-106 eV. 
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The values of the high energy cutoff of the incident spec- 
trum require some comments. Such values are clearly outside 
of the energy range of BeppoSAX instruments, but their de- 
termination is needed, in order to constrain accurately the re- 
flected spectrum. They seem not unreasonable. If the cutoff 
energy were inside the BeppoSAX energy range, the incident 
spectrum at higher energies would not provide enough photons 
to generate such an intense reflected component. The absence 
of a cutoff below 140 keV agrees with previous results (Dolan 
et al. 1 19871 Bouchet et al. 1993 1. In order to confirm the pres- 
ence of a high energy cutoff in the near MeV range further in- 
vestigations, e.g. with Integral, are required. 
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